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PROGRAM SUMMARY 



Title of the software package: First order ODE tools. 
Catalogue number: (supplied by Elsevier) 

Software obtainable from: CPC Program Library, Queen's University of Belfast, N. Ireland (see application form 
in this issue) 

Licensing provisions: none 

Operating systems under which the program has been tested: UNIX systems, Macintosh, DOS (AT 386, 486 and 
Pentium based) systems, DEC VMS, IBM CMS. 

Programming language used: Maple V Release 3/4. 

Memory required to execute with typical data: 8 Megabytes. 

No. of lines in distributed program, including On-Line Help, etc.: 4183. 

Keywords: 1°' order ordinary differential equations, symmetry methods, symbolic computation. 
Nature of mathematical problem 

Analytical solving of 1*^' order ordinary differential equations. 
Methods of solution 

Computer algebra implementation of Lie group symmetry methods. 
Restrictions concerning the complexity of the problem 

Besides the inherent restrictions of the method (there is as yet no general scheme for solving the associated PDE 
for the coefficients of the infinitesimal symmetry generator), the present implementation does not work with 
systems of ODEs nor with higher order ODEs. 

Typical running time 

This depends strongly on the ODE to be solved, usually taking from a few seconds to a few minutes. In the tests 
we ran (with 466 1=' order ODEs from Kamke's book |[, see seckh, the average times were: 6 sec. for a solved 
ODE and 20 sec. for an unsolved one, using a Pentium 133 with 64 Mb. RAM on a Windows 3.11 platform. 

Unusual features of the program 

The f^* order ODE-solver here presented is an implementation of all the steps of the symmetry method solving 
scheme; i.e., when successful it returns a closed solution, not the symmetry generator. Also, this solver permits 
the user to (optionally) participate in the solving process by giving an advice (HINT option) concerning the 
functional form for the coefficients of the infinitesimal symmetry generator (infinitesimals). All the intermediate 
steps of the symmetry method solving scheme are available as user-level commands too. For instance, using the 
package's commands, it is possible to obtain the infinitesimals, the related canonical coordinates, and the most 
general T^' order ODE invariant under a symmetry group. The package also includes a command for classifying 
ODEs (according to Kamke's' book|^) popping up Help pages with Kamke's advice for solving them, facilitating 
the study of a given ODE and the use of the package with pedagogical purposes. 
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LONG WRITE-UP 



Introduction 

The automatic computation of symbolic solutions for ODEs has already been implemented in almost all 
symbolic computation systems available. However, though Lie's discovery of symmetry methods (SM) put 
most of the known solving methods under a common umbrella[|l], ^, only a few of the computer algebra 
standard solvers use SM in their solving schemes and almost none use them for solving P* order ODEs. 

Generally speaking, one of the reasons for this is that the use of SM requires the solving of a system 
of partial differential equations (PDEs) for the coefhcients of the infinitesimal symmetry generators (in- 
finitesimals), whose solution is not obvious, thus involving the programming of semi-systematic heuristic 
methods for obtaining it^. Also, though in the specific case of 1^* order ODEs we just need to solve a 
single I''* order PDE, one may think that SM are useless anywayQ, since to guess a solution for that PDE 
would not be simpler than to guess an equivalent integrating factor for the original ODE. Furthermore, 
there are many well known methods, based on a previous classification of the ODE, which already give 
good results without using heuristics Q. So what would be the advantage of introducing heuristics in the 
solving process? 

There are however two strong arguments in favor of heuristic methods: they are the only chance of 
solving non- classifiable ODEs, and they are the main way to find the underlying logic structures with 
which one could develop new methods. 

Bearing all this in mind, we programmed an I''' order ODE-solver based on SM, implementing a set 
of heuristic methods for solving the PDE which leads to the infinitesimals. In addition, in order to test 
ODE-solvers in general, we prepared a file with the Maple input of the 576 ODEs| of Kamke's book, and 
another one with 424 ODES, mainly of Riccati, Abel and non-classifiable types. 

The set of user-level routines here presented also allows the user to obtain almost all the intermediate 
results involved in the SM scheme, to find the most general invariant ODE under a given symmetry, to 
test explicit/implicit results obtained by any (Maple) ODE-solver, and to classify ODEs, according to 
Kamke's book, including the optional pop up of Help pages with Kamke's advice about how to solve 
them. 

The goals of such a work can be summarized as: 

1. to build an ODE-solver complementing the standard dsolve, mainly for tackling non-classifiable 
ODEs, 

2. to build a research environment for discussing possible connections between patterns of infinitesimals 
and patterns of ODEs, mainly using the HINT option of the solver (see sec.^ in connection with 
the advisor (odeadv) command, 

3. to build a friendly educational environment for studying ODEs and their solving methods in general. 

The exposition is organized as follows. In sec.|l|, the SM scheme for solving ODEs is briefly reviewed 
and discussed. In sec.^ the package of routines is presented. Sec.^ contains a more detailed explanation 
of each of the heuristic methods implemented. Sec.^ is devoted to testing the package, mainly the solver, 
against the aforementioned set of 1,000 ODEs. These tests allow an evaluation of the solver from different 
angles, and include a comparison of the performances of the new solver and the Maple dsolve in solving 
Kamke's ODEs. Finally, the conclusions contain a brief discussion about this work and its possible 
extensions. 



1 Symmetry methods for 1 order ODEs 

Generally speaking, the key point of Lie's method for solving ODEs is that the knowledge of a (Lie) group 
of transformations which leaves a given ODE invariant may help in reducing the problem of finding its 
solution to quadratures |^. Despite the subtleties which arise when considering different cases, we can 
summarize the computational task of using SM for solving a given ODE, say, 

^From here on, we will use the abbreviate ODE for representing l**' order ODEs. 
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^ = Hy,x), (1) 

as the finding of the infinitesimals of a one-parameter Lie group which leaves Eq.(^ invariant, i.e., a pair 
of functions {^{y,x), 77(2/, x)} satisfying 

— + — -] ^ -L<|)2_^ „ — ^0 (2) 

dx \dy dx J dy dx dy 

followed by either 

a) the determination of the canonical coordinates, say {r, s}, of the associated Lie group, to be used in 

a change of variables which will reduce Eq.(^ to a quadrature; or, alternatively, 

b) the setting up of an integrating factor for Eq.(^ using the explicit expressions found for {S,{y, x), r](y, x)}. 

A first look at the symmetry scheme may lead to the conclusion that finding solutions to Eq. (||) would 
be much more difficult than solving the orig inal Eq.(0). However, what we are really looking for is a 
particular solution to Eq.(H), as simple as possible (for instance two constants), and for a non-classifiable 
ODE this particular solution may be the simpler one to be obtained. 

As an example, consider 

This equation is of Riccati type and cannot be solved easily^, but a polynomial guess for the infinites- 
imals (made by one of symgen's internal routines) leads to]^ 

> symgenC); # input = ODE, output = the infinitesimals 



J,{x,y)=x, jq{x,y)^x + y 

Passing Eq. (^) directly to odsolve (the solver) , it will call symgen, use the result above to determine 
the canonical coordinates, change the variables, and solve Eq.(^) as follows: 

> odsolve ("") ; 

xy/^ , ( ^/51n(a;) y/Z_Cl\ , , , x 
2^=__tanh( -^ + __ 1 +xln{x) + - 

What is amazing, and characteristic of symmetry methods is that changing {y — x ln(x))/ x to F{{y — 
x ln(x))/a;) in Eq.(||), where F is an arbitrary function, the symmetry remains unchanged and the solving 
scheme succeeds as well: 

^^ = F(y—^^\+Hx) (4) 
dx \ X I 



> odsolveC) ; 



y — X la{x) 



J (l + _a-F(_a)) 



The integral above is expressed using the new intat command (of the new version of PDEtools Q) 
which represents the integral evaluated at a point -analogous to the derivative evaluated at a point, intat 

^An analysis of Eq. (^) can be obtained by sending it in REDUCE format (to be tackled with the CONVODE program) 
to convode@riemann.physmath.fundp.ac.be 

*In what follows, the input can be recognized by the Maple prompt >. 
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displays the evaluation point as an upper limit of integration. For Eq.(^, Maple V R.4 returned an 
explicit result actually equivalent to the above (see Eq.(15)). 

Even when the built-in heuristics fail, in the framework of SM, it is possible to provide computational 
routines permitting users to test their own conjectures concerning the functional form of the infinitesimals 
(HINT option, see sec||). As an example of this, consider 



dx 



X + cos(e^ + {x+l) e"'') 



(5) 



This ODE is not classifiable nor is it solved by the automatic heuristics of odsolve. An example of 
a user-conjecture concerning the infinitesimals would be: 



^{x,y) = fix), 7j{x,y) = xg{y) 



(6) 



where f{x) and g(y) are unknown functions of their arguments. The HINT algorithm we programmed, 
given such an ansatz, will try to determine / and g, in this case leading to: 



> symgen(",HINT=[f (x) ,x*g(y)]) : 



> odsolve ("",HINT=[f(x) ,x*g(y)]) ; 



— e^, jq = xe 



y = In 2 arctan 



- {x + l)e~ 



2 The package's commands 

A detailed description of the package's commands, with examples and explanations concerning their 
calling sequences, is found in the On-Line Help. Therefore, we have restricted this section to a brief 
summary, followed by a detailed description of the solver, odsolve, and the advisor, odeadv, and a 
description of just one paragraph for each of the other routines^. Simple input/output examples can be 
seen in sec.||. 

2.1 Summary 

A brief review of the commands of the package is as follows: 





Command 


Purpose: 




odsolve 


solve ODEs using the symmetry method scheme 




fatint 


look for an integrating factor 




canoni 


look for a pair of canonical coordinates 




symgen 


look for a pair of infinitesimals 




equinv 


look for the most general ODE invariant under a given symmetry 




buisym 


look for the infinitesimals given the solution of an ODE 




odepde 


return the PDE for the infinitesimals 




odetest 


test explicit/implicit results obtained by any Maple ODE-solver 




symtest 


test a given symmetry 




odeadv 


classify ODEs and pop up a Help-page with Kamke's solving advice 



^This section contains some information already presented in the previous section; this was necessary to produce a 
complete description of the package. 
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2.2 Description 
Command name: odsolve 

Feature: V*' order ODE-solver based on symmetry methods 
Calling sequence: 

> odsolve (ode) ; } 

> odsolve(ode, y(x), way=xxx, HINT=[exprl, expr2] , int_scheme) ; } 

Parameters: 
ode - a 1"* order ODE. 

y(x) - the dependent variable (required when not obvious). 

way=xxx - optional, forces the use of only one (xxx) of the 6 internal algorithms {abacol, 

2, 3, 4, 5, abaco2} for determining the coefficients of the infinitesimal symmetry 
generator (infinitesimals) . 
HINT = [el,e2] - optional, el and e2 indicate a possible functional form for the infinitesimals. 
int_scheme - optional, one of: fat, can; meaning respectively: use an integrating factor; use 

canonical coordinates. 
Optional parameters can be given alone or in conjunction, and in any order. 
Synopsis: 

Given a 1*** order ODE, odsolve's main goal is to solve it in two steps: first, determine the infinitesi- 
mals of a 1-parameter symmetry group which leaves the ODE invariant, and then use these infinitesimals 
to integrate the ODE, by building an integrating factor or reducing the ODE to a quadrature using 
the canonical coordinates of that group. To determine the infinitesimals, odsolve makes calls to sym- 
gen, another command of the package, and to make the change of variables introducing the canonical 
coordinates, it makes calls to dchange, from the PDEtools package 

odsolve does not classify the ODE before tackling it, and is mainly concerned with non- classifiable 
ODEs, for which the standard dsolve fails. By default, odsolve starts off trying to isolate the derivative 
in the given ODE, then sequentially uses all six internal algorithms of symgen to try to determine the 
infinitesimals (see sec.^, and finally tries the two integration schemes sequentially. 

When odsolve succeeds in solving the ODE, it returns, in order of preference: an explicit closed solu- 
tion, an implicit closed solution, or a parametric solution (only when the derivative cannot be isolated^). 

All these defaults can be changed by the user. The main options she/he has are: 

1. to request the use of only one of the six internal algorithms for determining the infinitesimals 
(way=xxx option; different algorithms might lead to different symmetries for the same ODE, perhaps 
turning the integration step easier); 

2. to enforce the use of only one of the two alternative schemes for integrating a given ODE after finding 
the infinitesimals (can or fat optional arguments, useful to select the best integration strategy for 
each case); 

3. to indicate a possible functional form for the infinitesimals (HINT=xxx option). This option is 
valuable when the solver fails, or to study the possible connection between the algebraic pattern of 
the given ODE and that of the infinitesimals. 

4. to enforce the use of dsolve at first and symgen's algorithms only when it fails (assign the global 
variable sym := false); 

A brief description of how the HINT=xxx option can be used (see examples in sec.|l| and ||) is as follows: 

• HINT=[el,e2] indicates to the solver that it should take el and e2 as the infinitesimals and deter- 
mine the form of any (a maximum of two) indeterminate functions entering el and/or e2, such as 
to solve the problem (e.g., Eq.(^)). 

• HlNT=[el, '*'] where el represents the user's guess for the first infinitesimal, ^, and indicates 
to the solver that it must consider the other infinitesimal, 77, as a product of two indeterminate 
functions of one variable each: the independent and the dependent variables, respectively (e.g., 
Eq.([6l)). 
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• HINT=[el, ' + '] works as in the previous case, but replacing the product of two indeterminate 
functions by a sum of them . 

• HINT=parametric indicates to the solver that it should look for a parametric solution for the given 
ODE. This solving scheme may be useful even when the system succeeds in isolating the derivative 

0. 

Finally, there are three global variables managing the solving process, which are automatically set 
by internal routines but can also be set by the user, as desired. They are {dgun, ngun, sgun}, for 
setting, respectively, the maximum ctegree of some polynomials entering ansatze for the infinitesimals, 
the maximum number of subproblems into which the original ODE should be mapped, and the maximum 
size permitted for such subproblems. The dgun variable is automatically set each time symgen is called, 
according to the given ODE, whereas, by default, ngun and sgun have their values respectively fixed 
to 2 and 1 . Although the setting of these variables by the user might increase the efficacy of the 
algorithms significantly, increasing by one unit the ngun or sgun variables, depending on the given ODE, 
may geometrically slow down the solving process. 

Apart from the solver, the package includes an analyzer, odeadv, which can play a pedagogical 
role, or be a comfortable tool for classifying ODEs while studying the possible links between symmetry 
methods and standard ones. A detailed description of it is as follows: 

Command name: odeadv 

Feature: T** order ODE-analyzer and solving advisor 
Calling sequence: 

> odeadv(ode) ; 

> odeadv(ode, y(x) ,{typel,type2,type3, . . .},help) ; 

Parameters: 
ode - a V*- order ODE 

y(x) - the indeterminate function (necessary when not obvious) 

-[typel,type2, . . .> - optional, a subset of ODE classification types to be checked 
help - optional, to request the pop-up of a Help-page indicating Kamke's book 

advice for solving the given ODE. 

Synopsis: 

Given a 1"* order ODE, odeadv's main goal is to classify it according to Kamke's book and pop up a 
Help-page indicating Kamke's advice for solving it when required (by entering the word help as an extra 
argument). The Help-pages include examples and the Maple input lines implementing the advice, so as 
to allow the user to adapt them to her/his problem. 

When used without extra arguments, odeadv will try to classify the given ODE into one or more of 
the following types: 

quadrature; separable; linear: classes A, B, and C; homogeneous: classes A, C and D; exact; 
rational; Clairaut; Bernoulli; Riccati, Riccati special; Abel: 1"* type, 2^'^ type subclasses A, 
B and C; dAlembert; and "patterns": y = g{y'), x = g{y'), = G{x,y'), — G{y,y'), 
y = G{x,y'), X = G{y,y') 

The matching of the types is checked sequentially, in the order displayed above, odeadv may return 
more than one type when the ODE is of type homogeneous, exact or rational; otherwise, the first matching 
of a pattern will interrupt the process and a classification will be returned. 

As an option, the user can indicate the checking of just a subset of the types mentioned above, by 
giving this subset as an extra argument. 

A compact description of the purpose of the other routines of the package is as follows: 
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• symgen looks for the infinitesimals ^(x, y) and 77(2:, y) of a 1-parameter Lie group which leave the 
received ODE invariant, as a pair of functions satisfying Eq.(^). The options available for this 
command are the same way=xxx and HINT=xxx options of odsolve. 

• fatint and canoni respectively return the integrating factor for the given ODE or a set of trans- 
formations from the original coordinates to the canonical coordinates of the Lie group mentioned 
above. 

• equinv takes as arguments a list of two algebraic expressions, to be seen as the infinitesimals of a 
1-parameter Lie group, and returns, within the possibilities of the system, the most general ODE 
invariant under that group. 

• buisym takes as arguments the solution of an ODE and an indication of which are the dependent 
and independent variables, and looks for a pair of infinitesimals for the ODE which generated the 
problem. 

• odetest tests explicit and implicit solutions found by any Maple ODE-solver, such as dsolve or 
odsolve, by making a careful simplification of the ODE with respect to that solution. 

• symtest tests the results returned by symgen, returning when these results check ok, or an 
algebraic expression obtained after simplifying the PDE for the infinitesimals according to symgen's 
result. 



buisym and equinv may be useful in connection with the odeadv command and the HINT option of 
odsolve, in order to study the relation between symmetry patterns and ODE patterns. Also, it is worth 
mentioning that, since any solver will make extensive use of the whole Maple library, testing commands 
such as odetest and symtest would also be useful in detecting possible wrong results of the Maple 
library. 



3 The heuristic methods 
3.1 General remarks 

All the implemented heuristic methods assume that the system was successful in isolating the derivative, 
as in Eq.(^. These methods consist of six computational algorithms for building and testing an explicit 
algebraic expression for the infinitesimals and 77. By heuristic we mean that the possible success of the 
scheme cannot be determined a priori. 

Two of the six algorithms are in turn subdivided into eight and four subalgorithms respectively; thus, 
in all, the given ODE is actually tackled by trying sixteen different functional forms for ^ and 77. Three 
of these sixteen schemes involve different types of polynomial ansatze for 77}, whose form, apart from 
the polynomial degree, is virtually independent of the ODE input. In the other thirteen schemes, the 
ansatze for ^ and 77 depend on the received problem, and are not a priori of polynomial type, but rather 
are obtained by solving auxiliary ODEs. 

Four of the six algorithms directly tackle Eq. (||) , while the other two reformulate the problem taking 
into account that Eq. ^ admits 77 = ^ <I> as a general (useless) solution. Therefore, it is always possible 
to introduce 



^-e^ + X (7) 

in Eq.(^, where x(a;, y) is an unknown function of its arguments, in order to map the original problem 
into that of solving 

(8) 



dx dy \ dy 

for X- The knowledge of x, in turn, allows one to set ^ and 77 as desired using Eq.(0). 

For the case in which all these schemes fail, we programmed a HINT option to allow the user to try 



her/his own heuristics for finding the infinitesimals, as explained in sec. 2.2 
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The six algorithms 

The I''' algorithm, called abacoi, consists of 2 sets of 4 schemes each, leading to 8 different trials. The 
first 4 trials look for 77} taking one of the infinitesimals as and the other one as a function of only 
one variable, namely: 

U - 0, - /(x)}, {e = 0, 77 = /(y)}, {e = /(x),77 = 0}, ^ = /(y), ry = 0} (9) 

The possible success of each case is determined by algebraic factorization of Eq.(^), and when this 
possibility is confirmed, an auxiliary ODE is built to determine the explicit resulting form for /. 

The second sequence of 4 trials looks for 77} taking one of the infinitesimals as and the other 
one as an expression containing both x and y, namely: 



U = 0, r; = fix) g{y)}, = f{x) g{y), rj = 0}, 

{C^O,v^fiy)g{x)}, {e-/(y)g(x),77 = 0} (10) 

where g is an algebraic expression built by extracting factors from the numerator and denominator of 
the ODE to be solved, and / is an unknown function to be determined by solving auxiliary ODEs as in 
the first 4 trials. 
Examplf^: Kamke's ODE 120. 

> symgenC') ; 



_^ = 0, ji] = — e ^y (5(2/) = y s-nd f{x) was determined as — e ^) 

> odsolveC") ; 

y = e-n 

As a general benchmark, abacoi can determine a pair of infinitesimals for 151 of the 466 non- 
quadratures found in Kamke's first 500 ODEs, spending on each an average time of 0.26 sec. when 
successful and 1.1 sec. otherwise. 

The 2"'^of the set of six main algorithms consists of a bivariate polynomial in {x,y) ansatz for x to 
solve Eq.(H). The maximum degree of the polynomial can be set by the user by assigning the global dgun 
variable; otherwise, it will be determined automatically by an internal routine. This second algorithm 
can determine a pair of infinitesimals for 96 of Kamke's 466 non-quadrature ODEs mentioned above, 
spending an average time of 0.08 sec. when successful and 0.5 sec. otherwise. 

The 3'''^ and 4"^ of the set of six algorithms respectively consist of bivariate polynomial and rational 
ansatze in {x,y) for {^, ry}, to solve Eq.(||). 
Example: Kamke's ODE 236 (Abel, 2"'^type, class B), 

x{y + A) ^-y^-2y-2x^0 (12) 
ax 

> symgenC) ; 

_^ = 4 X + , jq = xy + lix 

In solving the aforementioned 466 Kamke's non-quadrature ODEs, these two algorithms can determine 
a pair of infinitesimals in 234 cases, spending an average time of 0.7 sec. in each success and 1.2 sec 
otherwise. 

The 5'^ algorithm builds an explicit expression for x{x,y), to solve Eq.(||), as follows: 

^In all the following examples, an alias has been set for y, via > alias(y=y(x)) ; 
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1. a basis of functions and algebraic objects is built by taking, from the given ODE, all the known func- 
tionsQ and composite algebraic objects, together with their derivatives, as well as all the unknown 
functions; 

2. a polynomial of degree 2 in such objects is built; its coefficients, in turn, are polynomials of degree 
d (the dgun variable mentioned above) in {x, y} with undetermined coefficients. 

This ansatz for x is introduced in Eq. (j^) resulting in a system of algebraic equations for the undeter- 
mined coefficients mentioned above. 
Example: Kamke's ODE 357 (no classification) 

^ ( ~ri ln(a;) sin(i/) + cos(y) (1 — a; cos(y)) ~ (13) 



> symgenC) ; 



dx 



4-0, 



cos(2y)a; x 



2 J 2x ln(x) sin(2/) 



> odsolveC") ; 



y = arccos 



2 ln(. 



X] 



2X + .C1 



This algorithm can find a pair of infinitesimals for 218 of Kamke's 466 non-quadratures mentioned 
above, spending on each an average time of 2.1 sec. when successful and 1.8 sec. otherwise. 

The 6**^ algorithm, called abaco2, leads to a sequence of 4 trials for the infinitesimals 77}, making 
use of the HINT option of symgen, with two indeterminate functions of only one variable each, namely: 

{£. - 9{x), V = f{x)}, {C = g{x), r] = /(y)}, 

{^^9{y),v^f{x)}, {^^9{y),v = .f{y)} (14) 

The possible success of these trials cannot be determined by algebraic factorization of Eq.(^), as in 
abacoi, except in a few cases. Hence, the alternative scheme we implemented for managing such HINTs 
consists of determining / and g as follows: 

1. subdivide Eq.(^ into subexpressions involving only one of {/, g}; 

2. build a list of candidates for / and for g with the solutions to these subexpressions; 

3. build a list of pairs of candidates by taking one candidate from each list. 

To avoid spending much time, the number of candidates of each list in step (^) is restricted, by default, 
to n = 2, leading to a maximum of 4 pairs of candidates for each functional form in Eq. (p^) . This default 
can be changed by the user by assigning the global variable ngun. 
Examples: ODE 593 (file 6) with an arbitrary function F{x, y) 

dx y/y \ ' 

> symgenC) ; 



4=4' (giy) ^ fix) ^ ^ - seeEq.^) 

e"" y/y ^ 

> odsolveC'"); # Maple V R.4 output using intat and RootOf 

^By known function we mean a function whose derivative rule is known by the Maple system. 
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y = ( RootOf ( j 




> odetestC'"" , ") ; # explicit and implicit results can be tested with odetest 



For this ODE, Maple V R.3 returned an implicit (equivalent) result, as in Eq.(Q). 

abaco2 can determine rj} for 114 of Kamke's 466 non-quadratures mentioned above, spending on 
each an average time of 0.4 sec. when successful and 8 sec. otherwise. 

To conclude this section, a curious result illustrating the HINT option is given by the infinitesimals 
found by symgen for the general form of Bernoulli's equation, 

^ + f{x)y + g{x)y'^ = 0, (16) 
where / and g are arbitrary functions of their arguments, and a is a symbolic power. The giving of a 



HINT, asking for a separation of variables by product for one of the infinitesimals (see subsec.2.2) leads 
to: 

> symgenC ,H1NT=[0, '*'] ) ; 



4 = 0, _77 = e 
> factor (odsolve ( " " , H1NT= [0 ,'*'])) ; 



f{x)dx (a~l) 



f{x)dx {a — 1) 
y = \ e^-^ / (a — 1) 



f{-a)d-a (a — 1) 
.g(_a)e \ J / d.a ~ _C1 



4 Tests 



We tested the set of routines here presented having two different ideas in mind: to confirm the correctness 
of the returned results, and to evaluate the new routines' performance in solving Kamke's ODEs. The 
performance test was done with the intention of comparing the efficacy of a SM-based solver such as 
odsolve with that of a classification-based solver as is Maple's dsolve. 

We used 10 test files, containing a set of 1,000 ODEs distributed in: files 1 to 6 containing the 576 
ODEs of Kamke's book, and files 7 to 10 containing 424 ODEs of Riccati, Abel and non- classifiable 
types, for which the classification-based dsolve (R.3/4) fails. The ODEs of files 7-10 were built using the 
equinv command, departing from symmetries of Kamke's ODEs, and were used to test the correctness 
of the results returned by symgen and the integration subroutines of odsolve. The "fail by dsolve" 
condition for these 424 ODEs was proposed to evaluate odsolve in solving the type of ODEs for which 
it was really designed, i.e., non-classifiable or usually not solved by dsolve. 

The purpose of preparing files with such a big set of ODEs was to set up what would be a convenient 
test-file for evaluating the performance of computational implementations in solving classifiable/ non- 
classifiable ODEs. The 10 files, with the input code in Maple, Reduce and Mathematica format, are 



available at our site: http://dft.if.ucrj.br/symbcomp.html 



A summary with the number of ODEs of each file belonging to each classification (we used our odeadv 
classification command) is as follows: 
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Class: 


File: 


Totals: 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


quadrature 


10 





1 


12 


10 


13 














46 


separable 


19 


10 


3 


8 


10 


7 














57 


linear 


14 


15 


16 


1 


1 

















47 


homogeneous 


1 


19 


15 


11 


28 


4 


14 


19 


19 


8 


138 


exact 








14 


14 




















28 


rational 


5 


48 


68 


28 


26 


1 


17 


64 


54 


40 


351 


Clairaut 











7 


12 


8 














27 


Bernoulli 





1 


























1 


Riccati 


23 


46 


1 











13 


27 


19 


21 


150 


Abel 


15 


8 


30 











12 


19 


25 


42 


151 


d'Alembert 


3 


2 


15 


7 


36 


10 


1 











74 


none 


16 


11 


6 


27 


23 


36 


91 


27 


25 


42 


304 


has F{x, y) 


21 


4 


8 


8 


1 


13 


43 








7 


105 


Total of ODEs 


100 


100 


100 


100 


100 


76 


124 


100 


100 


100 


1,000 



Table 1. Classification of the 1,000 ODEs of files 1-10 



The quadratures pointed out in the table were found after isolating the derivative, receiving a rhs 
not containing x or y{x) (the independent and dependent variables all along the 1,000 ODEs). The 
next-to-last line of the table above indicates the number of ODEs having an arbitrary function F{x,y); 
e.g. F(x'^ — /y). Also, in many cases, odeadv returned more than one classification for a given ODE; 
all these possible classifications are reflected above. 

4.1 Test of symgen 

As explained in sec.^ the SM approach for solving ODEs involves two main steps: the determination 
of a pair of infinitesimals and their subsequent use in the integration process. We therefore divided the 
tests into the same two steps. 

The first performance test, thus, concerned the explicit determination, by the symgen command, of 
infinitesimals for Kamke's ODEs (only the 466 non-quadratures of the first 500 ODEs were used). The 
table below displays the total number of successes, the average time spent with each solved/fail case, and 
the number of successes of each of symgen's six algorithms when the whole set of ODEs was tackled 
using only one of them. The second column of the table indicates the number of non-quadrature ODEs 
per file. 



File 


ODEs 


Successes 


Average time^ 


5 and rj can be determined by 


ok 


fail 


abacoi 


2 


3 


4 


5 


abaco2 


1 


90 


44 


1.7 sec. 


4.2 sec. 


33 


13 


21 


25 


29 


23 


2 


100 


71 


1.2 sec. 


2.0 sec 


30 


13 


47 


46 


39 


27 


3 


99 


81 


0.7 sec. 


3.6 sec. 


35 


36 


55 


46 


40 


15 


4 


88 


72 


3.8 sec. 


4.1 sec. 


38 


28 


44 


45 


49 


21 


5 


89 


75 


1.7 sec. 


5.1 sec. 


15 


6 


67 


71 


61 


28 


Totals: 


466 


343 


Si 8 sec. 


« 12 sec. 


151 


96 


234 


233 


218 


114 



Table 2. Kamke's ODEs for which the infinitesimals were determined by symgen; 73% 



4.2 Test of odsolve's integration schemes and comparison of performances 

The two integration alternatives we programmed for the second step of the solving scheme consist of either 
using an integration factor or making a change of variables (canonical coordinates) reducing the ODE to 
a quadrature. Of course, the corresponding routines work only if symgen succeeds in determining the 
infinitesimals. 

Also, as mentioned, though odsolve was thought as a complement to the Maple dsolve, we prepared a 
comparative test of performances to have an idea of the possible efficacy of a SM-based solver, if compared 
to that of a classification-based solver as dsolve. With this purpose, we tested the standard dsolve with 
the non-quadrature Kamke's ODEs too. The results we found in testing the integration schemes and in 
comparing odsolve's performance with that of dsolve can be summarized as follows: 

*Due to the large number of ODEs being tested, we interrupted and computed as fail all calculations consuming more 
than a few minutes. This was the case in Kamke's ODEs numbers 51, 340, 347, 373, 394, 460 and 479; they were recalculated 
separately, receiving an answer in 5 ^ 20 min. for all but ODE 479 . 
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File 


ODEs 


Solved by 


Successful integration via 


dsolve 


odsolve /sy mgen 


can 


fat 


1 


90 


47 


44/44 


41 


40 


2 


100 


71 


71/71 


66 


70 


3 


99 


75 


79/81 


61 


78 


4 


88 


63 


70/72 


51 


68 


5 


89 


69 


70/75 


69 


55 


Totals: 


466 


325 


334 / 343 


288 


311 



Table 4. Comparative performance 



It was a surprise to us that, even when more than 80 % of the ODEs used in the test are of classifiable 
type (see Table 1.), the performance of odsolve, which does not classify the received ODE before tackling 
it, was almost the same as that of dsolve^. Also, odsolve's integration schemes succeeded in 97% (334 
of 343) of the cases in which symgen determined a pair of infinitesimals, convincing us in that the main 
problem is, actually, the determination of those infinitesimals, not their posterior use in the integration 
process. 



5 Conclusions 

This paper presented a computer algebra implementation of symmetry methods for solving 1^' order 
ODEs. Despite the additional complications introduced by the necessity of solving an auxiliary 1^' order 
PDE, this method proved to be a valuable tool for solving non-classifiable (or even classifiable) ODEs, 
as shown in sec.^ resulting in the extension to Maple's dsolve we were looking for. 

Moreover, we would like to remark the interactive character of the solver too; i.e., the user is given 



the possibility of participating in the solving process (the HINT option, see sec. 2. 2), achieving in this 
way a significant extension of the solving capabilities of the scheme. 

Also, the computational routines presented here were designed not only for solving ODEs: the odeadv 
command and the HINT option of odsolve were thought as a combination for investigating new solving 
methods, and the set of commands here presented was intended to be complete from the educational point 
of view. Actually, with them it is possible to tackle an ODE in the framework of symmetry methods, 
directly or step by step, it being also possible to go back using the equinv and buisym commands, 
completing the conceptual cycle. 

On the other hand, this is a first version of the package and many things can be improved. To 
mention but a few, the classification abilities of odeadv can be extended, and used by symgen in order 
to determine the most efficient ordering of the algorithms for tackling the given ODE (at present this 
order is unmovable); also, a study of the fail examples can be used as a basis for building new algorithms, 
mainly for ODEs of Riccati and Abel types. Such new algorithms can be implemented straightforwardly 
using the HINT option of symgen. Moreover, for instance, it is possible to extend the ideas presented 
here to tackle higher order ODEs too. These subjects are part of our present interests and we expect to 
be able to report related work in the near future 0. 
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